Internal Josephson-Like Tunneling in Two-Component Bose-Einstein Condensates 
Affected by Sign of the Atomic Interaction and External Trapping Potential 
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We study the Josephson-like tunneling in two-component Bose-Einstein condensates coupled with 
microwave field in respond to various attractive and repulsive atomic interaction under the various 
aspect ratio of trapping potential and the gravitational field. It is very interesting to find that the 
dynamic of Josephson-like tunneling can be controlled from fast damped oscillations and asymmetric 
occupation to nondamped oscillation and symmetric occupation. 
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The existence of a Josephson current is a direct man- 
ifestation of macroscopic quantum phase coherence [lj 
and have numerous important applications in condensed 
matter physics, quantum optics and cold atom physics, 
for example, precision measurement, quantum computa- 
tion. The physical origin of the Josephson current is the 
temporal interference of the two systems which must both 
have a well defined quantum phase and a different aver- 
age energy per particle, respectively. Recently, the ex- 
perimental realization of multi-component Bose-Einstein 
condensates (BECs) of weakly interacting alkali atoms 
has provided a route to study Josephson effect @, 0, 0, 0, 
IS 0)11, EH in a controlled and tunable way by means of 
Feshbach resonance [IJ [H, [H, Q [H [H, [jj, M, [H, 
so far unattainable in charged systems. The physical 
origin of the low-energy Feshbach resonances is the low- 
energy binary collisions described by the difference of the 
initial and intermediate state energies which can be ef- 
fectively altered through variations of the strength of an 
external magnetic field. 

In this Letter, we study the dynamics of Josephson-like 
tunneling in two-component Bose-Einstein condensates 
coupled with microwave field under tuning the atomic in- 
teraction from attractive case to repulsive case in the var- 
ious aspect ratio of trapping potential. Through chang- 
ing the relative value between attractive intra-atomic and 
inter- atomic interaction, we find that, occupation on rela- 
tive number of atoms represents two very different behav- 
ior, which one is fast damped oscillation, but the other is 
nondamped and irregular oscillation. In particular, when 
the atomic interaction is changed from attractive to re- 
pulsive case, we find that relative number of atoms also 
represents nondamped and irregular oscillation. 

We consider a system of longitudinal elongated two- 
component BECs coupled by the microwave field with 
the effective Rabi frequency Q and finite detune 8, where 
the Q and 8 are independent upon time and space coordi- 
nate as experimental case Q. At zero temperature, this 
system can be described by macroscopic wave function ipj 
(normalized to unity, i.e., f(\ipi\ 2 + \ip2\ 2 )dz = I) which 
obey the dimensionless one-dimensional Gross-Pitaevskii 
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FIG. I: Temporal evolution of relative number of atoms 
Nd for initial Gaussian wave function in response to vari- 
ous aspect ratio of trapping potential 7. (a) Correspond- 
ing to attractive atomic interactions for the parameters as 
Gn = G22 = 4Gi2 = -10 and fi = l.5u±, 8 = 0.5lo ± . 
(b) Corresponding to attractive atomic interactions for the 
parameters as Gn = G22 = 4Gi2 = —10 and Q = 1.5u)±, 
8 = -0.5c.jx. 



equations (GPEs) as [1, 0] 
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where j, k = 1, 2, j 7^ k, and w±_ is radial harmonic trap- 
ping frequency, time and coordinate unit are 2/lo± and 
a± = \J h/rnuj±, respectively, and V(z) is the longitudi- 
nal confining potential, Gjj — (4NHajj)/(maj_u)±) and 
Gjk = (ANhajk) j '(ma^_w±) are dimensionless intra- and 
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FIG. 2: Temporal evolution of relative number of atoms 
Nd for initial Gaussian wave function in response to vari- 
ous aspect ratio of trapping potential 7. (a) Correspond- 
ing to attractive atomic interactions for the parameters as 
4G U = 4G 22 = G12 = -10 and Q = 1.5w±, 5 = 0.5lu ± . 
(b) Corresponding to attractive atomic interactions for the 
parameters as 4Gn = AG22 = G12 = —10 and fl = 1.5u>x, 
8 = -0.5o> x - 



FIG. 3: Temporal evolution of relative number of atoms 
Nd for initial Gaussian wave function in response to vari- 
ous aspect ratio of trapping potential 7. (a) Correspond- 
ing to repulsive atomic interactions for the parameters as 
Gn = G22 = 4Gi2 = 10 and Q, = 1.5wj_, 5 = -0.5wx. 
(b) Corresponding to repulsive atomic interactions for the 
parameters as 4Gn = 4G22 = G12 = 10 and Q, — 1.5u)±, 
8 = -0.5lo±. 



inter-species effective atomic interaction potentials, re- 
spectively, where djj and ajk are corresponding s-wave 
scattering lengths, respectively, N is the total atomic 
number, and m is the atomic mass assumed the same 
for both two species. 

In the following discussion, we consider two-hyperfine 
spin states of 7 Li bose atom which the magnitude and 
sign of the atomic interaction can be tuned to any value 
via a magnetic-field dependent Fcshbach resonance. And 
we choose dimensionless V(z) = j 2 z 2 + Xz, where 7 = 
^z/^_l is aspect ratio of trapping potential and A deter- 
mines the strength of the gravitational field, lj z is longi- 
tudinal harmonic trapping frequency and ui± = 2tt x 625 
Hz (so, the time and coordinate units are 0.51 ms and 
a±_ = 1.51 /im respectively). Moreover, considering 
the condensates will be unstable and collapse when the 
number of particles is large enough (l7| . so we choose 
N = 6000, which provides a safe range of parameters for 
avoiding instability and collapse occurring. 

First of all, we consider the attractive case, i.e., intra- 
and inter-species effective atomic interaction potentials 
are time-independent attractive case, which can be easily 
realized using Feshbach resonance experimentally. And 
the effective Rabi frequency and finite detune are cho- 
sen as Q = 1.5o;j_ and 8 = 0.5cjx corresponding to Fig. 
la and i5 = — 0.5wx corresponding to Fig. lb. Here, 
we have chosen an = 022 is just considering the experi- 



mental feasibility and assumed that tuning atomic inter- 
action wouldn't destroy the condensates in this system. 
We denote relative number of atoms between two species 
as Nd, i-e., N d = Ni - N 2 , where ATj = ffy^dz and 
^2 = J\ip2\ 2 dz, respectively. We study Josephson-like 
tunneling in this case and the results are shown in Fig. 
1, which is based on numerical simulation of Eqs. ([T]) 
for the initial Gaussian wave function under the various 
aspect ratio of trapping potential. As shown in Fig. 1, 
during the evolution process, exchange of atoms is always 
damped very fast and almost suppressed completely irre- 
spective of aspect ratio of trapping potential. In particu- 
lar, it is important to notice that occupation on the each 
components represents very asymmetric behavior and the 
atoms almost transform to second component completely 
in the case of 8 = 0.5o;_l, but almost transform to first 
component completely in the case of 8 = — 0.5lu±. In this 
case, system has lower energy through exchanging the 
atoms from first to second components under the con- 
dition 8 = 0.buj±, but from second to first components 
under the condition 8 = — 0.5uj±. 

From the experimental points of view, it is very im- 
portant to study the dynamics of Josephson-like tunnel- 
ing under the tuning atomic interaction which provides a 
route to study Josephson effect in a controlled and tun- 
able way by means of Feshbach resonance so far unattain- 
able in charged systems. So, we consider the attractive 
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FIG. 4: Under the gravitational field, temporal evolution of 
relative number of atoms Nd for initial Gaussian wave func- 
tion in response to various aspect ratio of trapping potential 
7. (a) Corresponding to attractive atomic interactions for 
the parameters as Gn = G22 = 4Gi2 = — 10 and Q — 1.5lo±, 
S = 0.5uj±. (b) Corresponding to attractive atomic interac- 
tions for the parameters as 4Gn = AG22 = G12 = — 10 and 
= 1.5uj±, S — Q.5u±. 



intra- and inter-species effective atomic interaction po- 
tentials become 4Gn = AG22 = G12 — — 10. The results 
are shown in Fig. 2, which is based on numerical simu- 
lation of Eqs. fl]) for the initial Gaussian wave function 
under the various aspect ratio of trapping potential. We 
can see that amplitude of oscillation becomes larger with 
increasing the aspect ratio of trapping potential 7. Es- 
pecially, contrary to first case, occupation on the each 
components almost represents very symmetric behavior 
and the atoms are not favorite to stay on only one of two 
components because the energy of system in asymmetric 
occupation is larger than symmetric occupation on the 
each components. 

Next, we study the dynamics of Josephson-like tunnel- 
ing under the repulsive atomic interaction in response to 
various aspect ratio of trapping potential. As shown in 
Fig. 3(a) corresponding to the case Gn = G22 = 4Gi2 = 
10 and Q = 1.5uj±, S = —0.5lj±, amplitude of oscillation 
becomes larger with increasing the aspect ratio of trap- 
ping potential 7, and occupation on the each components 
also represents symmetric behavior approximately, which 
is like the case shown in Fig. 2. However, oscillation be- 
comes slower with increasing the aspect ratio of trapping 
potential 7, which is very different from the case shown 
in Fig. 2. It is very interesting to notice that Fig. 3(b) 
corresponding to the case 4Gn = 4G22 = G12 = 10 rep- 
resents very different phenomena from Fig. 3(a). 




FIG. 5: Under the gravitational field, temporal evolution of 
relative number of atoms Nd for initial Gaussian wave func- 
tion in response to various aspect ratio of trapping potential 
7. (a) Corresponding to attractive atomic interactions for 
the parameters as G11 = G22 = 4Gi2 = 10 and fl = 1.5u)±, 
S = 0.5lo±. (b) Corresponding to attractive atomic interac- 
tions for the parameters as 4Gn = 4G22 = G12 = 10 and 
Q = 1.5lj±, S — 0.5ui±. 



Another interesting case is to study the dynamics of 
Josephson-like tunneling under the gravitational field, so 
V(z) in Eqs. |T]) becomes 7 2 z 2 + Xz, where A determines 
the strength of the gravitational field and equals 0.839 
in unit of (l/2)hui±. As shown in Fig. 4(a) correspond- 
ing to attractive atomic interaction and Gn = G22 = 
4Gi2 = —10, = 1.5lu±, S = 0.5lu±, we can see that the 
relative number of atoms is damped slower than the case 
without the gravitational field when 7 is very small, but 
is damped with the same speed as the case without the 
gravitational field when 7 is large as shown in Fig. 1. In 
Fig. 4(b) corresponding to attractive atomic interaction 
and 4G U = 4G 22 = Gi 2 = -10, n = 1.5w_l, S = 0.5uj ±1 
we can see the same phenomena as the case without the 
gravitational field as shown in Fig. 2, except the ampli- 
tude and frequency of oscillation of relative number of 
atoms have some differences. 

When the atomic interaction becomes repulsive under 
the gravitational field, the dynamics of Josephson-like 
tunneling is shown in Fig. 5. We can see that under the 
atomic interaction satisfied Gn = G22 = 4Gi2 = — 10 
corresponding to Fig. 5(a), the dynamics of Josephson- 
like tunneling under the gravitational field represents the 
same phenomena as the case without the gravitational 
field as shown in Fig. 3(a), except the amplitude and 
frequency of oscillation of relative number of atoms have 
some differences. But under the atomic interaction sat- 
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isfied 4Gn = AG22 = G12 = 10 corresponding to Fig. 
5(b) and 7 is very small, the dynamics of Josephson-like 
tunneling under the gravitational field is very different 
from the case without the gravitational field as shown in 
Fig. 3(b). 

In conclusion, we find a robust way to control the 
dynamics of Josephson-like tunneling from fast damped 
oscillations and asymmetric occupation to nondampcd 
oscillation and symmetric occupation in two-component 
Bose-Einstein condensates coupled with microwave field. 
Recent developments of controlling the scattering length 
and realization of multi-component BECs in the exper- 
iments allow for the experimental investigation of our 
prediction in future. 
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